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Abstract 


The Badhwar- O’ Neill (BON) Galactic Cosmic Ray (GCR) model is based 
on GCR measurements from particle detectors. The model has mainly 
been used by NASA to certify microelectronic systems and the analysis 
of radiation health risks to astronauts in space missions. The BON14 
model numerically solves the Fokker-Planck differential equation to ac- 
count for particle transport in the heliosphere due to diffusion, convec- 
tion, and adiabatic deceleration under the assumption of a spherically 
symmetric heliosphere. The model also incorporates an empirical time 
delay function to account for the lag of the solar activity to reach the 
boundary of the heliosphere. This technical paper describes the most 
recent improvements in parameter fits to the BON model (BON14). Us- 
ing a comprehensive measurement database, it is shown that BON14 is 
significantly improved over the previous version, BON11. 
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1 Introduction 


It is well documented in the literature that the ionizing nature of Galac- 
tic Cosmic Ray (GCR) particles poses a potential health risk for crew 
members in space, particularly for future long-term missions in free- 
space [1-3]. Long-term exposure to the mixed (both low and high) Linear 
Energy Transfer (LET) GCR radiation significantly increases the Risk of 
Exposure Induced Cancer (REIC) and Risk of Exposure Induced Death 
(REID) [1-3]. As with exposure to any charged particle, GCR particles 
may interact with the cell at the deoxyribonucleic acid (DNA) level and 
cause DNA damage, and free-radical particle production. 

Another major concern arises from the interaction of the GCR par- 
ticles with the electronics inside and outside of a spacecraft. Energetic 
GCR particles may deposit energy in electronics, e.g., microprocessors, 
memory units, sufficient to cause memory bit flips and latch-up, which 
are called Single Event Effects (SEE) [4,5]. 

In order to assess astronaut exposure and electronic effects associ- 
ated with GCR exposure, models are used to describe relevant spectra 
of particles and energies appearing in deep space. These models formu- 
late the GCR transport, associated energy loss mechanisms, which are 
also known as GCR modulation (either interstellar or in the solar sys- 
tem), and flux as a function of energy spectrum. The Badhwar- O’ Neill 
(BON) model [6-9] numerically solves the Fokker-Planck (FP) equation 
that takes into account diffusion, convection, and adiabatic deceleration 
within the heliosphere. The solution to the FP equation as implemented 
in the BON model is the flux of GCR particles of a given charge, Z, 
as a function of energy. The solution is obtained under the assump- 
tions of a quasi-steady state and a spherically symmetric interplanetary 
medium [10,11]. Under these assumptions, the FP equation can be writ- 
ten as [12]: 


r 2 dr 


(r 2 V s lI) 




]_d_ 

j . 2 



) (i) 


where r is the radial position in astronomical units (AU); T is the kinetic 
energy (MeV/n); U(r,T) is the GCR flux (particles/sr-m 2 -s-MeV/n), 
V’s(r) the solar wind speed (~ 400 km/s); n(r,T) the particle diffusion 
coefficient tensor; and a(T ) = (T + 2Tq)/(T + To), with To being the 
rest energy of the GCR particle. The solution also assumes that at 
a boundary distance r = Rb , modulation of U(r,T) is negligible, and 
therefore provides the boundary condition at U(Rb,T ) = To as a known 
quantity. This quantity, Uq, is ion specific and parametrically described 
by several free parameters, which are known as Local Interstellar (LIS) 
parameters. 

The analytical solution of the full FP equation exists only for a simple 
form of k, and with constant values of V s (r) and a(T) [13]. Therefore, 
the FP equation must be solved numerically to obtain flux in the desired 
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energy range and solar conditions [12]. The solution depends on the 
assumption that several key quantities are known, including the set of 
parameters known as Local Inter-stellar (LIS) parameters. 


2 Initial parameters and data sets 

2.1 Local Interstellar parameters 

As described earlier, at a distance well outside of the solar system, the 
GCR modulation due to the turbulent solar wind and heliospheric mag- 
netic field is negligible. At around Ri = 100 AU, it is assumed that each 
GCR ion energy-flux spectrum is constant. This constant GCR field is 
referred to as the local interstellar spectrum ( ji on ) and represents one 
of the boundary conditions for the BON model. The LIS equation is 
described with the following parametric model: 

Uo oc Jion(T) U 6=i00 ^= JO P S (T + To) -7 (2) 

where jo, 6, and 7 are free parameters for each particular GCR ion, and 
(3 = v/c is the velocity of the ion relative to the speed of the light. The 
LIS parameters are formulated by using the GCR measurement data 
from detectors at or near 1 AU, e.g., satellite and balloon measurements. 
In the model, the flux of any ion beyond nickel (Z > 28) is obtained by 
scaling from the silicon result. 

In this paper, we describe the most recent LIS parameters fitted to an 
updated data set to improve the BON model. It should be emphasized 
that all GCR energy spectra (measurements and model) reported in this 
paper apply to free-space beyond the Earth’s magnetosphere. 

2.2 Description of the LIS parameters in BON14 model 

In this revision of the BON model, we have modified the LIS parameters, 
jo, 5, and 7, driven by a sensitivity study by using several metrics. For a 
detailed description of the sensitivity analysis, please see the referenced 
publication series by T. Slaba et al. [14 16]. This analysis was completed 
by dividing ions into charge and energy groups as shown in Table 1. 

One of the most distinguishing results of this study was that for the 
differential effective dose rate as a function of kinetic energy behind 20 
g/crn 2 of aluminum shielding at a period of minimal sun activity, the 
two energy groups between 0.5 GeV/n and 4.0 GeV/n account for most 
of the exposure for all ions. As seen in the table, for instance, GCR 
protons and alphas with boundary energy less than 0.5 GeV/n induce 
approximately 9% of the total effective dose. GCR ions with Z > 2 
and boundary energies less than 0.5 GeV/n induce less than 4% of the 
total effective dose behind 20 g/crn 2 aluminium shielding. Based on the 
results of this study, the new LIS parameters are fitted to the GCR data. 
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The measurement database and uncertainty metrics described in T. 
Slaba et al. [14 16] were implemented and automated. This allowed 
a range of parameter combinations to be evaluated through the BON 
model to determine which parameters lead to minimal uncertainties be- 
tween the model and measurement database for each energy and ion 
group. A similar procedure, to calibrate free parameters in a simplified 
empirical GCR model, was previously described and utilized by Matthia 
et al. [17]. 

In the BON case, the uncertainty was quantified using average rela- 
tive difference between data and model, relative uncertainty on effective 
dose, and y 2 analysis. Readers are referred to Eqs.[l-6] in T. Slaba et 
al. [16] for a detailed description of how the relative uncertainty results 
were incorporated as weighting factors in the analysis. 

The computation of the average uncertainty in each energy bin al- 
lowed either positive or negative uncertainty values in each bin, which 
permitted to systematic model trends to be seen (negative values corre- 
spond to underprediction, and positive values correspond to overpredic- 
tion) . Whereas the final parameter set for each ion was chosen somewhat 
subjectively, reduced uncertainty results in the upper, most relevant en- 
ergy bins were preferred. 


Table 1. Relative contributions of the indicated charge groups to effective 
dose behind 20 g/cm 2 of aluminum shielding at a period of minimal 
solar activity. The ions were grouped into 5 energy and atomic number 
categories. Charge groups: Zi =1, Z2 = 2, Z3 = 3 — 10, Z4 = 11 — 20, and 
Z5 = 21 — 28. Kinetic energy groups: Ei= 0.0 — 0.25, E2 = 0.25 — 0.5, 
E3 = 0.5 — 1.5, E4 = 1.5 — 4.0, and E5 = 4.0 — 1.0 x 10 4 GeV/n. 


Energy 

Charge 

Ei 

E2 

E3 

e 4 

e 5 

Total [%] 

Zi=l 

1.2 

5.4 

18.2 

18.4 

14.8 

58.1 

Z2=2 

1.2 

2.2 

4.1 

2.9 

1.7 

12.2 

Z 3 =3 - 10 

<0.1 

3.3 

3.8 

1.3 

0.8 

9.1 

Z 4 =ll - 20 

<0.1 

0.2 

6.6 

2.0 

1.1 

10.0 

Z 5 =21 - 28 

<0.1 

<0.1 

4.7 

3.8 

2.1 

10.6 

Totals 

2.5 

11.1 

37.4 

28.4 

20.5 

100.0 


The LIS parameters for each of the corresponding ions from H to Ni 
are shown in Table 2. 
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Table 2. The LIS parameters that are used in the B0N14 model shown 
for each ion. 


z 

7 

4' 

jo 

1 

2.75 

-2.82 

9.50x 10~ 4 

2 

2.80 

-2.00 

4.53x 10 -5 

3 

3.21 

-0.69 

6.37xl0 -8 

4 

2.93 

1.50 

1.20x 10 -7 

5 

3.00 

-0.40 

2.40x 10 -7 

6 

2.70 

-2.00 

1.60x 10 -6 

7 

2.95 

-0.60 

2.65x 10 -7 

8 

2.73 

-1.90 

1.50x 10 -6 

9 

3.08 

0.40 

1.63x 10 -8 

10 

2.75 

-1.60 

2.35x 10 -7 

11 

2.73 

-1.80 

4.60x 10 -8 

12 

2.70 

-2.40 

3.03x 10 -7 

13 

2.75 

-1.40 

5.30xl0 -8 

14 

2.65 

-2.40 

2.65xl0 -7 

15 

3.15 

2.00 

5.68xl0 -9 

16 

2.70 

-1.00 

5.78x 10 -8 

17 

3.13 

2.00 

5.99x 10 -9 

18 

2.90 

0.60 

1.68x 10 -8 

19 

3.13 

0.80 

7.90x 10 -9 

20 

2.75 

-1.60 

3.23x 10 -8 

21 

3.15 

0.40 

3.50x 10 -9 

22 

3.00 

-0.50 

1.44xl0 -8 

23 

3.00 

-0.50 

7.14xl0 -9 

24 

2.90 

-1.00 

1.78xl0 -8 

25 

2.80 

-1.00 

1.39x 10 -8 

26 

2.60 

-2.40 

2.00x 10 -7 

27 

2.60 

-2.50 

l.llxlO -9 

28 

2.55 

-2.40 

1.19x 10 -8 
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2.3 Selection of GCR data 

In the past, the BON model parameters were uniquely influenced by 
measurements from the Cosmic Ray Isotope Spectrometer (CRIS) on the 
NASA Advanced Composition Explorer (ACE) spacecraft. The CRIS 
instrument is currently measuring the flux of ions and their isotopes 
from boron (Z=5) to nickel (Z=28). 

Although the lowest and highest kinetic energy measurements are 
ion specific, the CRIS provides kinetic energy of GCR isotopes between 
~ 50 - 500 MeV/n overall. In the BON14 model, greater emphasis was 
placed on the higher kinetic energies, a region not covered by ACE/CRIS. 
Nonetheless, we show that the updated model still accurately represents 
the lower energy regions. 

In this revision of the BON model, BON14, we have included the 
GCR data beyond 1970, which spans Solar Cycle 20 to 24 (to date). 
Table 3 lists the name, flight type, measurement times, range of GCR 
ion charge and energy of the measurement data sets that we used for 
calibrating the LIS parameters. 


Table 3. The list of GCR data sets that were used to calibrate the 
LIS parameters in the BON model. A similar table with more detailed 
information can be found in [16]. 


Name 

Flight 

Time 

Ions (Z) 

Energy [GeV / n] 

ACE/CRIS [18] 

Satellite 

1998 - Present 

5 - 28 

0.05 - 0.5 

AMS [19,20] 

STS-91 

1998 

L 2 

0.1 - 200.0 

ATIC-2 [21] 

Balloon 

2002 

1 - 26 f 

4.6 - 10.0 

BESS [22] 

Balloon 

1997 - 2000, 2002 

1, 2 

0.2 - 22.0 

CAPRICE [23, 24] 

Balloon 

1994, 1998 

L 2 

0.15 - 350.0 

CREAM-II [25] 

Balloon 

2005 

6 - 26+ 

18.0 - 10.0 

HEAO-3 [26] 

Satellite 

1979 

4 - 28 

0.62 - 35.0 

IMAX [27] 

Balloon 

1992 

1, 2 

0.18 - 208.0 

IMP-8 [28] 

Satellite 

1974 

6 - 14+ 

0.05 - 1.0 

LEAP [29] 

Balloon 

1987 

L 2 

0.18 - 80.0 

MASS [30] 

Balloon 

1991 

L 2 

1.6 - 100.0 

PAMELA [31,32] 

Satellite 

2006 - 2009 

1, 2 

0.08 - 10.0 

TRACER [33] 

Balloon 

2003 

8 - 26+ 

0.8 - 10.0 

Lezniak et al. [34] 

Balloon 

1974 

4 - 26 f 

0.35 - 52.0 

Minagawa et al. [35] 

Balloon 

1975 

26, 28 

1.3 - 10.0 

Muller et al. [36] 

STS-51 

1985 

6 - 14+ 

50.0 - 10.0 

Simon et al. [37] 

Balloon 

1976 

5 - 8 

2.5 - 10.0 


' Not all ion data are available in the provided range. 


Most of the GCR data that are collected by the ACE/CRIS instru- 
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ment have very low count rates, which leads to high statistical mea- 
surement errors. There are a few abundant ions with low statistical 
measurement uncertainty for ions, such as B, C, N, O, and Fe. However, 
almost all of the ions have very high statistical uncertainty, such as data 
provided for F, P, Cl, Sc, and most of the heavy ions. Hence, the tables 
that report the comparison of BON model metrics, multiple successive 
Bartels rotations were averaged to reduce the statistical uncertainty in 
the measurements. By doing this, most of the ions have an uncertainty 
of less than 10% and the highest uncertainty, in this case for Co, was 
reduced to 25%. 

3 Solar modulation 

One of the important input parameters of the BON model is the decelera- 
tion potential, </>, (a.k.a. solar modulation parameter). The deceleration 
potential is correlated with the diffusion coefficient, k, and may be used 
to account for attenuation of the LIS due to the state of the heliosphere. 
Larger values of <p represent lower GCR flux. The <f> is in units of rigid- 
ity (MV), and is proportional to the momentum/charge required for a 
particle to penetrate the heliosphere. The parametric relation between 
k and 4> can be described by the following equation [38]: 

«(»■, t) = ^/LP[ 1 + ( — ) 2 }/Ht) (3) 

v s r 0 

where V s is the constant solar wind speed (~ 400 km/s), r is distance 
from the Sun in AU, t is time in years, kq is a diffusion constant (1.6x 10 21 
cm 2 /s), and P is particle rigidity in units of MV. In order to fit ions from 
H to Ni by using a simple form of the LIS model, ?’o is arbitrarily set 
to 4.0 AU [38]. It is important to note that the value of /to is chosen 
somewhat arbitrarily, therefore the values of </>(f) are not significant. 
They simply define the level of modulation for a given period of time 
within the context of our choice of parameters. In Fig. 1, GCR data for 
various ions for low- and high-solar activity times are shown. The solar 
modulation is very noticeable for energies lower than ~ 4 GeV/n for H 
and He. 

Methods that determine (j) from the current measurement of solar 
activity at the Sun, such as sunspot number, show a lag of the GCR 
modulation. The BON model incorporates a time-delay function [39] to 
account for the lag because of Solar effects which take a certain amount 
of time (~ 8 — 14 months) to reach the boundary of the heliosphere. In 
addition, the magnetic polarity reversal of the Sun is considered by the 
model. 

On the other hand, direct measurement of the GCR flux by a space- 
craft instrument samples the current state of the heliosphere. The sunspot 
method has the advantage of predicting future GCR fluxes; however, the 
spacecraft method can more directly sample the actual GCR flux. 
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Energy [MeV/n] 


Figure 1. The flux data as a function of kinetic energy for H, He, 0, 
and Fe are shown for all the data used for calibration (this plot excludes 
ACE/CRIS data). 


In Fig. 2, the <j> computed by using BON14 and monthly smoothed 
average international sunspot number (SSN) [40] as a function of time 
are shown. The left axis represents the SSN and the right axis represents 
the 4>. The figure clearly presents the time-delay effect incorporated in 
the cj), where 4> is calculated from the SSN by accounting for the time- 
delay between the time of the solar activity (the sunspots) and the time 
this effect (the magnetic field disturbance) reaches into the heliosphere 
enough to significantly modulate the GCR flux. 

In this current revision of the BON model, as opposed to the previous 
BON11 model [9] where there were two cj) values (one for the proton and 
another 0 for all other ions); there is only one < j> input for all of the ions. 
A single value of 4> describes the state of the heliosphere at a given time 
and is used for all elements from hydrogen (Z=l) through plutonium 
(Z=94). 

Analysis shows that the correlation of the spectra of all the significant 
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Figure 2. The monthly smoothed average SSN and (deceleration po- 
tential) as a function of time for the period between 1998-2015 are shown. 
On the left axis SSN (red solid line) and on the right axis <f> (blue dashed 
line) are shown. The cj) function is shifted in time due to the time-delay 
function. 

GCR elements (Z=l to 28) is better when using direct sampling of the 
heliosphere by spacecraft to determine the 0. This is important to the 
SEE analyst who needs to know the actual history of GCR fluxes, where 
the SEE analyst needs knowledge of actual worst-case conditions for 
spacecraft designs. 

It is worth noting here that the <j) is directly associated with solar 
activity and can be determined (in principle) by a number of methods. 
The various solar indices - e.g., sunspot number, radio and X-ray flux, 
neutron monitor rates, spacecraft GCR measurements, and even tree 
growth rings - indicate the temporal variation of the Sun’s intensity. 

4 Comparison of the BON14 model to the data 

The comparison of the BON14 model with the GCR measurement data 
(GCR database provided in Table 1) was done by using several metrics 
as explained in Sec. 2.2. That included looking at relative differences, 
analyzing distributions of uncertainties, and inclusion of measurement 
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0[MV] 


uncertainty. It is worth mentioning here that since almost all of the 
GCR measurement are broadly binned (with error bars), we integrated 
the BON model output data to match the same energy bin of the mea- 
surement data and considered the average of the points. No significant 
differences were found if comparisons were done using the data mid- 
points instead. 

The relative difference, Rd, is a measurement between the GCR mea- 
surement data and the BON model that allows us to find where the model 
underpredicts or overpredicts the measured GCR flux. Rd is shown as: 


„ , 1 Modelk — Datak 

Rd = — > — 

N Datak 

k = 1 


( 4 ) 


where N is the number of measurements. The average of the residual 
Rd does not provide a complete picture of the model versus measure- 
ment agreement, as the quantity may end up in the vicinity of zero due 
to positive-negative cancellation effect, yet the uncertainty may still be 
high. Therefore, we also used the average absolute relative difference, 
\Rd\, to determine the overall difference between the model and the GCR 
measurement data, which can be described as: 


q ” \Model k - Data,] 

11 N Data k V ' 

This metric better quantifies the spread in model errors and may be 
used as a quality check for the model. In addition to the nominal data 
measurement comparison metrics, we extended our uncertainty analy- 
sis by including a measure of fit relative to the error bars described as 
following: 


_ 1 -A Modelk - [Data k - £ k ] 

U lower — AT / v - 5 W 

N Dat(lk ~ £ k 

TT _ 1 Modelk - [Data k + £ k ] 

•upper Datak+£ + ’ 

where and e + are published lower and upper error bars on data, re- 
spectively, and Ui ower and U upper represent the uncertainty interval that 
incorporates the GCR data measurement error bars. The width of the 
interval is approximately considered as the measurement uncertainty. In 
Fig. 3, the comparison of uncertainty metrics between the GCR measure- 
ment database presented in Table 3 and BONIO [8], BON11, and BON14 
models are shown as functions of energy and charge groups. As clearly 
seen from the figure, there is a noticeable improvement with BON14 in 
almost every bin over BONIO and BON11. The improvement is espe- 
cially significant for ions in the important energy bins between 0.5 and 
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4.0 GeV/n. As shown in Table 1, based on the sensitivity study, about 
66% of the total effective dose is induced by the ions within this energy 
range. Thus, more emphasis was placed on minimizing the uncertainty 
in those bins. 



Energy (GeV/n) 


u 


upper 




<=|Rd 


^ lower 


Figure 3. The average relative model uncertainty as a function of charge 
and energy groups are shown for BONIO, BON11, and BON14 models. 
As illustrated on the right side of the figure, the marker positions repre- 
sent Rd, while error bars represent upper and lower uncertainty interval. 


Figures 4—7 show the comparison of the average of the integrated 
differential flux for H, He, O, and Fe for each data set with BON11 and 
BON14 models. In these figures, we separated the data into two energy 
regions to simplify comparisons. The top plots represent lower energy 
bins (< 4.0 GeV/n) and bottom plots represent higher energy bins (> 
4.0 GeV/n). We define the average flux (< U >) as the integral of 
the differential flux ( dU/dE ) over all energy bins, which is given by the 
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following equation: 


<U>= 


1 

Nbin 


dU Irn 
dE iE ' 


( 7 ) 


where iV&j n is the number of energy bins and dE is the bin width of 
the data set. Average flux was evaluated by summing over the reported 
energy bins. The \Rd\ values for low-energy for H ion are 18.7% and 
13.6% for B0N11 and BON14 models, respectively. Moreover, as shown 
in the Fig. 5 for He, Fig. 6 for O, and Fig. 7 for Fe, the BON14 model 
presents a significant improvement over the BON11 model especially for 
energies below 4.0 GeV/n. 

In Fig. 8 Rd comparisons between ACE/CRIS data with the BON11 
model (top) and BON14 model (bottom) were plotted for O [at a rec- 
ommended] mid-energy of 230.8 MeV/n. The left axis on this figure 
presents the GCR measurement data and BON model differential fluxes 
in units of particles/s-sr-m 2 -MeV/n. The right axis shows average Rd 
for each data set. The \Rd\ for BON14 is about a factor of 2 lower than 
BON11; where it was calculated as 26% for BON11 and 13.6% for B014 
for the oxygen GCR data between 1998 - 2014 (August through Decem- 
ber, 1997 data were not included in the validation due to issues with 
data calibration as reported in [41]). 

As discussed in Sec. 2, greater emphasis was placed on higher kinetic 
energies, which is not covered by the ACE/CRIS model. Although the 
new model still has large relative difference during solar maximum pe- 
riods (~ 2000 — 2005), we present that BON14 has been significantly 
improved over the BON11 version. Moreover, BON14 follows the oxygen 
GCR data better for the solar activity period between 2010 and 2015. 
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Figure 4. The average flux for H ions as a function of measurement 
date is shown for BON11 and BON14 models for low-energy (top) and 
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ACE/CRIS. 
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Figure 8. The ACE/CRIS flux data as a function of Bartels’-rotation (27- 
day average) from 1998-present for 0 ion is compared with (top) BON11 
model and (bottom) BON14 model. The axis on the right represents Rd 
between the BON model and GCR data. 
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In Table 4, the overall average relative difference metrics, Rd and 
\Rd\, are presented for the BON11 and BON14 models. These values 
were determined by averaging the relative difference metrics over each 
ion, as given in Appendix A. The results in Appendix A were obtained 
by considering the entire GCR measurement database (all ion and en- 
ergy bins). It is seen in Table 4 (and Appendix A), that for BON11, 
the \Rd\ value is calculated as 23.7% and Rd = 17.9%. For BON14, 
these errors are reduced to \Rd\ = 13.0% and Rd = —0.4%. The BON11 
model systematically overpredicts the GCR measurement data, whereas 
the new BON14 model provides a more balanced prediction. More im- 
portantly, the overall spread in the new model error was reduced from 
23.7% to 13.0%. It should be emphasized that at all energies, BON14 
is only a marginally improved fit to the hydrogen GCR data compared 
with BON11. There is room for improvement since BON14 underpre- 
dicts the hydrogen GCR data during ~ 8 years since 1998 (per lower 
panel of Fig. 4). A similar underpredicting pattern is observed for He as 
well. 

Table 4. The average relative difference metrics between the entire GCR 
database with BON11 and BON14. 


All GCR measurement data 


Average [%] 


Rd 

\Rd\ 

BON11 

17.9 

23.7 

BON14 

-0.4 

13.0 


5 Conclusions 

Presented are the new methodology to obtain the LIS parameters and 
comparison of the results of the new BON14 model with a new set of GCR 
data from various balloon, satellite, and space shuttle measurements. 
The sensitivity study described earlier showed that the two energy groups 
between 0.5 GeV/n and 4.0 GeV/n account for most of the shielded 
effective dose for all ions. Therefore, we modified the LIS parameters, 
jo, 5, and 7 , driven by this sensitivity study. 

The overall average Rd value is calculated as 17.9% for the BON11 
and -0.4% for the BON14 model, whereas \Rd\ is found as 23.7% and 
13.0% for BON11 and BON14 models, respectively. The uncertainty 
metrics for 5 energy bins and charge groups also revealed that overall 
the BON14 model has been improved significantly with respect to the 
previous models both at low- and high-energy regions. The new model 
especially provides less uncertainty for ions in the important energy bins 
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between 0.5 and 4.0 GeV/n, which account for a significant part of the 
total effective dose behind shielding. However, it should be noted that 
the LIS parameters used by BON14 slightly underpredict hydrogen GCR 
data for most of the years since ~ 1998 (beginning of a solar maximum). 
Following the release of the BON14 model, we will attempt to improve 
the model by incorporating several physics effects on the GCR particles 
(e.g., curvature and gradient drifts) which have not been included in the 
previous BON models. 
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Appendix A 


The average relative difference metrics 


Table A5. The comparison of the metrics between the entire GCR 
database with BQN11 and BON14 for each ion. 


BON11 

BON14 


ACE/CRIS 

Others 

ACE/CRIS 

Others 

Z 

Rd 

\Rd\ 

Rd 

\Rd\ 

Rd 

\Rd\ 

Rd 

\Rd\ 

1 

- 

- 

7.5 

13.6 

- 

- 

-1.6 

11.1 

2 

Q 

- 

- 

18.3 

19.1 

- 

- 

0.3 

11.7 

O 

4 

_ 

_ 

61.5 

61.5 

_ 

- 

-7.3 

20.1 

5 

22.3 

24.4 

46.0 

46.0 

1.3 

11.1 

-1.5 

16.1 

6 

21.3 

23.6 

27.1 

34.2 

-2.1 

12.7 

0.0 

15.6 

7 

21.4 

23.3 

32.8 

35.0 

-5.0 

12.7 

1.4 

35.8 

8 

17.0 

20.6 

28.6 

34.4 

-3.1 

12.8 

-0.3 

13.7 

9 

14.3 

18.5 

24.3 

24.3 

8.6 

13.0 

-12.4 

14.8 

10 

15.1 

17.9 

16.4 

30.9 

2.9 

11.4 

-3.3 

13.5 

11 

13.8 

17.9 

40.9 

40.9 

3.1 

12.1 

-7.5 

24.1 

12 

16.6 

18.7 

22.9 

31.7 

5.5 

12.6 

-1.7 

17.1 

13 

14.2 

17.3 

33.8 

33.8 

6.9 

12.5 

-3.6 

16.9 

14 

14.4 

17.2 

26.0 

29.4 

2.9 

12.2 

-2.0 

16.8 

15 

6.6 

16.2 

19.8 

19.8 

7.8 

15.1 

-4.3 

4.3 

16 

11.0 

14.9 

15.1 

24.7 

3.3 

11.8 

-5.0 

12.2 

17 

2.5 

14.2 

22.5 

22.5 

3.9 

11.7 

-1.3 

3.0 

18 

5.7 

13.6 

-3.0 

18.5 

3.2 

11.9 

-6.3 

12.4 

19 

-0.6 

13.6 

16.1 

16.1 

10.0 

14.0 

-2.5 

2.5 

20 

-0.5 

14.2 

20.2 

20.2 

-9.5 

15.2 

-0.1 

15.3 

21 

-5.9 

16.0 

17.1 

17.1 

6.6 

13.7 

2.3 

4.3 

22 

4.5 

13.0 

25.5 

25.5 

-1.9 

11.7 

1.7 

3.9 

23 

2.2 

13.8 

23.3 

23.3 

5.5 

12.2 

2.6 

2.9 

24 

5.6 

13.5 

19.9 

19.9 

6.9 

12.9 

-2.5 

2.5 

25 

8.7 

15.1 

33.5 

33.5 

-0.1 

12.9 

-2.1 

2.5 

26 

7.1 

12.1 

28.5 

30.9 

-1.9 

12.5 

-1.6 

11.5 

27 

-11.1 

36.8 

2.9 

15.7 

10.1 

34.0 

-14.3 

14.3 

28 

2.7 

14.2 

44.4 

44.4 

-2.7 

13.4 

-0.4 

12.2 

Average 

8.7 

17.5 

24.9 

28.4 

2.6 

13.6 

-2.7 

12.3 
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Table A6. The comparison of the metrics between entire GCR data 
(ACE/CRIS and Others) with BON11 for each ion. The comparison pre- 
sented here is separated into low-energy (< 4.0 GeV/n) and high-energy 
(> 4.0 GeV/n) regions. The average relative difference is averaged over 
the results of all the ions. 


BON11 


ACE/CRIS Others All GCR Data 

Low-energy Low-energy High-energy All energies 


z 

Rd 

\ Rd \ 

Rd 

\ Rd \ 

Rd 

\ Rd \ 

Rd 

\ Rd \ 

1 

- 

- 

8.6 

18.7 

6.5 

8.4 

7.5 

13.6 

2 

q 

- 

- 

21.8 

22.0 

14.9 

16.2 

18.3 

19.1 

o 

4 

_ 

_ 

103.9 

103.9 

19.1 

19.1 

61.5 

61.5 

5 

22.3 

24.4 

63.4 

63.4 

28.5 

28.5 

34.1 

35.2 

6 

21.3 

23.6 

44.4 

44.4 

9.8 

24.1 

24.2 

28.9 

7 

21.4 

23.3 

37.9 

37.9 

27.7 

32.1 

27.1 

29.1 

8 

17.0 

20.6 

45.5 

45.5 

11.6 

23.3 

22.8 

27.5 

9 

14.3 

18.5 

31.5 

31.5 

17.0 

17.0 

19.3 

21.4 

10 

15.1 

17.9 

42.8 

42.8 

- 10.0 

19.1 

15.7 

24.4 

11 

13.8 

17.9 

67.9 

67.9 

14.0 

14.0 

27.3 

29.4 

12 

16.6 

18.7 

46.7 

46.7 

- 0.8 

16.6 

19.7 

25.2 

13 

14.2 

17.3 

49.9 

49.9 

17.6 

17.6 

24.0 

25.5 

14 

14.4 

17.2 

49.7 

49.7 

2.3 

9.2 

20.2 

23.3 

15 

6.6 

16.2 

19.7 

19.7 

19.9 

19.9 

13.2 

18.0 

16 

11.0 

14.9 

30.1 

30.1 

0.2 

19.3 

13.0 

19.8 

17 

2.5 

14.2 

25.9 

25.9 

19.1 

19.1 

12.5 

18.3 

18 

5.7 

13.6 

0.4 

14.0 

- 6.4 

23.0 

1.3 

16.0 

19 

- 0.6 

13.6 

19.4 

19.4 

12.9 

12.9 

7.7 

14.8 

20 

- 0.5 

14.2 

30.3 

30.3 

10.1 

10.1 

9.8 

17.2 

21 

- 5.9 

16.0 

15.0 

15.0 

19.2 

19.2 

5.6 

16.5 

22 

4.5 

13.0 

32.4 

32.4 

18.5 

18.5 

15.0 

19.2 

23 

2.2 

13.8 

35.5 

35.5 

11.2 

11.2 

12.8 

18.6 

24 

5.6 

13.5 

28.5 

28.5 

11.3 

11.3 

12.7 

16.7 

25 

8.7 

15.1 

47.6 

47.6 

19.4 

19.4 

21.1 

24.3 

26 

7.1 

12.1 

44.8 

44.8 

12.2 

16.9 

17.8 

21.5 

27 

- 11.1 

36.8 

18.6 

18.6 

- 12.9 

12.9 

- 4.1 

26.3 

28 

2.7 

14.2 

61.7 

61.7 

27.0 

27.0 

23.6 

29.3 

Average 

8.7 

17.5 

37.9 

38.8 

11.8 

18.0 

17.9 

23.7 


24 





Table A7. The comparison of the metrics between entire GCR data 
(ACE/CRIS and Others) with BON14 for each ion. The comparison pre- 
sented here is separated into low-energy (< 4.0 GeV/n) and high-energy 
(> 4.0 GeV/n) regions. The average relative difference is averaged over 
the results of all the ions. 


BON14 

ACE/CRIS Others All GCR Data 


z 

Low-energy 
Rd \Rd\ 

Low-energy 
Rd \Rd\ 

High-energy 
Rd \Rd\ 

All energies 
Rd \Rd\ 

1 

- 

- 

1.8 

13.6 

-5.0 

8.7 

-1.6 

11.1 

2 

O 

- 

- 

1.4 

12.9 

-0.8 

10.4 

0.3 

11.7 

O 

4 

_ 


2.2 

23.5 

-16.7 

16.7 

-7.3 

20.1 

5 

1.3 

11.1 

-4.0 

18.9 

1.1 

13.3 

-0.1 

13.6 

6 

-2.1 

12.7 

-6.8 
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